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Abstract. We present a method to semiclassically compute the pair creation rate of bosons and fermions in de 
Sitter spacetime. The results in the bosonic case agree with the ones in the literature. We find that for the constant 
electric field the fermionic and bosonic pair creation rate are the same. This analogy of bosons and fermions in the 
semiclassical limit is known from several flat spacetime examples. 


1 Introduction 

Quantum field theory (QFT) in curved spacetime is one way of merging Einstein’s theory of gravitation 
and QFT in the usual Minkowski spacetime within a self consistent framework |[T] [2) . One of its central 
results is the discovery of particle creation in a time-dependent gravitational field 0. This mechanism 
is believed to generate the primordial cosmic inhomogeneities that serve as seeds for the observed large- 
scale structure of the universe. 

A lot of studies of QFT in curved spacetime are investigating de Sitter (dS) spacetime. It has constant 
scalar curvature and is maximally symmetric for any given dimension, in the same way as Minkowski 
spacetime. In cosmology, dS space is used as a model for both the early stage of inflation (for reviews see 
||4l ) and the late stage of acceleration of the expansion (an introductory review is given in 0). In addition 
to its relevance for cosmology, dS space may give hints to understand the quantum nature of spacetime 
0. These facts motivate the investigation of physical processes in this spacetime. 

Originally, particle creation was studied under a strong field background originating from a time depen¬ 
dent vector potential 0. Since then, it has been explored in a more general set-up, where the effects of 
the gravitational field and the electric field are both contributing to the creation of pairs. The Schwinger 
mechanism in curved spacetime, more specifically, in dS space, has recently started to be studied in depth 
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Better understanding Quantum electrodynamics (QED) in dS spacetime could be insightful for a couple 
of problems. It is an interesting framework for the study of false vacuum decay or bubble nucleation 
mm. Constrains on magnetogenesis scenarios were put in 0 via the backreaction of the created pairs 
and its induced current. In some models of preheating, it might give clues on the open problem of baryo- 
genesis. Via the AdS-CFT correspondence, it has also been used as a playing field to test the ER=EPR 
conjecture Col . Additionally it might help to better understand renormalization schemes in curved space- 
time mm and the relations of these schemes to each other. 

The purpose of this paper is to investigate spin 1/2 pair production under the influence of an external 
electric field, in four-dimensional dS spacetime (dS 4 ). This problem has already been studied for bosons 


in two-dimensional dS spacetime (clS2 J ID and in dS4 (8). In lfT3lfl9l spin 1/2 pair production in dS2 was 
investigated but the 4 D analog case is much less studied. We propose a derivation of the number of pairs 
created in the semiclassical limit following the ideas of ED. We find that there is no difference between 
bosons and fermions in this limit. 

The paper is organized as follows: in section[ 2 ]we investigate the Schwinger effect in dS4. After describ¬ 
ing the basic equations, we show how the Dirac equation becomes a coupled second order differential 
equation in the presence of both electric and gravitational terms. We review the derivation of the bosonic 
pair creation rate and calculate the fermionic one with the help of a semiclassical expansion for general 
electric fields in dS spacetime in section [ 3 ] We use this result in section [ 4 ] to compute the bosonic and 
fermionic pair creation rate of a constant electric field in dS spacetime. Finally we draw some conclusions 
in section[ 5 ] 


2 Dirac field in dS Space-time 


We consider QED coupled to a Dirac fermion in dS4. In order to study Schwinger pair production, we 
assume that the gravitational and electrical field are background fields and that the fermionic field is 
dynamic. The action is given by 


S = 


f d 4 x y/-g£ = f d 4 xV=g 


K 2 


— mif/if/ - 




where the field strength is defined in the usual way F MV = d M A v - d v A fl . 
The dS4 spacetime we want to study is described by the metric 


(1) 


ds 2 = a(rj) 2 (dr ] 2 - dx 2 ), 


( 2 ) 


with signature Here 77 is the conformal time which is parametrized by the Hubble factor in the 

following way 
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( 3 ) 


For the description of spinors the tetrad field is used. It is related to the metric through the relation 


guv — e^e^ab. 


( 4 ) 


where r/ a b is the usual Minkowski metric. Throughout this paper, we use Greek indices for spacetime 
indices (//, v - rj, x, y, z) and Latin indices for tetrad ones (77, b = 77, x, y, z). Applying the tetrad formalism 
to dS4, one gets 
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( 5 ) 


The covariant derivative for fermion fields is defined as 


Vfj = hBfj + ieAfj(x) - cr ab . 


( 6 ) 


with the commutator of the gamma matrices 


Cab = 


'a, 7b] 
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and the spin connection are defined as 


of = ^ [e ba (x)d ll e a a (x) - e aa (x)d^ b a (x) + e aa (x)d a e b (x) - e ba (x)d a e alJ (x ) 
e bv (x)e aA (x)e c/1 (x)d A e c v (x) - e av (x)e bA (x)e cfl (x)d A e c v (x)]. 


The non-zero components of the spin connection ® can be shown to be 
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where prime denotes derivative with respect to conformal time rj. 

The gamma matrices V' are related to the gamma matrices in the tangent flat space y" 

yf_ = 7n e f 

We will work with the Dirac representation of the gamma matrices, i.e. 
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Varying the action with respect to the spinor field gives the Dirac equation 

(*y% - mj 77) = 0. 

Using Eqs. (0) and (O, this equation becomes 

|i {y 1 hd^ + —aHy° + ieA^Y^ - tj ) = 0. 


( 8 ) 
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One now considers the auxiliary field Ttx, rj) = a A ^ 2 (ri)i//(x, //) which can be thought of as the equivalent 
Mukhanov-Sasaki variable in inflation models JT|. With this substitution the Dirac equation takes 

rm 


of the 
the form 


{y^iihd/j - eAf,) - if T(x, 77 ) = 0 . 
We will also decompose this field in momentum modes according to 


^(x, 77 ) ~ e'' k V k . 


(14) 


(15) 


To solve the Dirac equation, for the purpose of the calculation of the pair creation rate, it is often useful 
to use the squared version of the Dirac equation because of its similarities to the Klein-Gordon equation, 
see e.g. mm. The squared Dirac equation can be found using 


T'tx, 77 ) = / \iidfj - eA^Ti)) + ma ( 77 )] 0 (x, 77 ) 
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in the Dirac equation. In the previous equation, 0k(?7) is the Fourier transform (in the sense of (fl~5l >) of 
(f>(x, rj). Here we consider a background vector potential for the electromagnetic sector such that 


= A(jj)5y 


( 18 ) 



For such fields the squared Dirac equation takes the form 


(19) 

( 20 ) 


( 'h 2 d 2 + u>]l(ti) 2 - ihma'( t])^ (p* ± iheA'{rf)<p 2 {if) = 0 , 

( h 2 d 2 + tj^(rj) 2 + ihma'(T]fj(f >2 ± iKeA'{rj)<p*{rj) = 0, 

where the effective pulsation and the kinetical momentum are defined as 

Uk(r]) 2 = p z (rj) 2 +k 2 L +m 2 a(r]) 2 , p z (r]) = k z + eA(rj), k\ = kr x + k 2 . (21) 

This equation can be compared to the equivalent bosonic problem. The equation of motion derived from 
the Klein-Gordon equation is (see Eq (2.13) of ED) 
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(i);/ + CJli hoson)//) j Pk — 0, 


tUk.boson( ? 7)~ — tUk(//)“ ^ • 

T 


( 22 ) 


(23) 


The equation of the bosonic problem (l20l can be understood as a harmonic oscillator with a time depen¬ 
dent pulsation. The two other terms in Eq. (l20l) are new and due to the fermionic nature of the particles 
considered. On the one hand, the mass term was already derived e.g. in 03 where no background electric 
field was considered. On the other hand, the electric term is present for instance in Il22l in flat spacetime. 
Such that this equation is a generalization of the Dirac equation in curved spacetime with background 
electric and gravitational field. 

The squared Dirac equation is also analogous to the Dirac equation for two-component fields in flat space- 
time. In ED) a method was used to semiclassically compute the pair creation rate for these fields. It is 
possible to use the same method for the case studied here. Instead of looking for a solution of the squared 
Dirac-equation (l20l we will however use the ansatz 
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This ansatz can be derived by finding the solution of the squared equation analogous to 1201 and then use 
to construct a solution for the Dirac equation (fT4t . Observe that ^^( 77 ) and ^^( 77 ) are independent 
since 


<Ak,t( 7 7) t • <Ak,i(f7) = 0 . (25) 

Putting (l24l > in the Dirac equation (ITD) leads to the equations 

ih ^'*( 77 ) + ma(Tj) (Af ( 77 ) ± (p z (r/) + i k ± ) ^( 77 ) = 0 , (26) 

ih ^(p) - ma{if) ij/tip) + (p-(p) - i k ± ) ^( 77 ) = 0, (27) 

which we will solve in section [3] 

3 Semiclassical number of pairs in dS spacetime 

The semiclassical expansion is considering cases where a vacuum state for the produced particles exits 
in the asymptotic future. That is true if the background fields are evolving slowly. To quantify the slow 
varying background more precisely, we introduce a dimensionless slowness parameter T by replacing 






the scale factor a(rf) by a family of functions «■/ (//) = a(jj/T). Doing that, in the limit of infinitely slow 
varying backgrounds, T —» oo, the derivatives of a(rj) will tend to zero. Orders of T are usually called 
adiabatic orders ID but it can be noticed that the only place where T is involved is in the derivative d n of 
the Dirac equation so it is possible to "formally" pose T — l/h and expand in power of h. 

In this section, we will first review the calculation of the number of particles in the bosonic case, which 
can be carried out with flat spacetime techniques and then present the fermionic case. The strategy will 
be the following: 

• Reformulate the equation of motion in terms of an equation for the mode functions a(rf), POf). 

• Perform a multiple integral iteration to compute |/?| 2 . 

• Calculate the integrals with a semiclassical saddle point approximation to derive the number of 
created pairs for each momentum mode k. 


3.1 Equations for the mode functions 

In this section we derive the equations for the mode functions from the Klein-Gordon and Dirac equation 
respectively. The aim is to construct equations in which a' (if) depends only on and vice versa, in 
order to perform the multiple integral iteration of the next section. 


3.1.1 Bosonic case 

To compute the semiclassical pair creation rate we start from (l22t . The form of these equations is the same 
as in flat spacetime, the only difference being the specific time dependence of the fields. We will shortly 
review the well know techniques of flat spacetime (see e.g. ifZTIO applied to our field configurations. One 
can use the ansatz (which is inspired by a WKB expansion) 
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where 


KqOt) - ^ J' cj k (r)dT, (30) 


where a(r]),/3(rj) are the mode functions. 

The momentum spectrum of the pair creation rate is defined as 

n k = lirn \J3(rj)\ 2 , (31) 

i}—> 0 


with the boundary conditions 


P (-°°) = 0, 


a(- 00 ) = 1. 


(32) 


Using this, together with the ansatz (I28l >- (l29t in the Klein-Gordon equation (l22l> we find that the mode 
functions are connected through coupled differential equations 
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3.1.2 Fermionic case 


To derive analogous equations in the fermionic case we start from the equations (l26l) and (ITH) . The initial 
idea of this work comes from the similarities between this equations and Eq. (11-12) of l20l . For the 
semiclassical treatment one makes the following ansatz 

tfW - - = : (a ± (7?)[ Wk ( ?7 ) - rnam** m +0*0 7)[«*fa) + > (35) 

V w k(?7) SPziri) -1 k± 

= — 7 = - r— + ma(q)]e~^ K(r,) + 78 * 07 ) [ 0 ^( 77 ) - ma{q)]^ m) \ , (36) 
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with the integrals 
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(38) 


Using the ansatz (l35l>- (l36l > in (l26l and (f27l) . we find that the mode functions are connected through coupled 
differential equations 
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representing the fermionic corrections to the analog bosonic case. 


3.2 Multiple integral iteration 

In this section we will perform the multiple integral iteration for the fermionic case. Because of the similar 
form of the equations (l33l - (l34l and (l39l- (l40l) the results for the scalar case can be derived from the ones 
obtained below by setting G a (rf) — Gp(rf) = 1 and Kit]) —> Kq(t]). 

By iteratively using Eqs. (l39l > and (l40t and the boundary conditions (l32t one finds 
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As described in |[23l these integrals are dominated by the classical turning points given by 

oJkiPp) = 0 . 


(44) 




















It can be shown (see H2QI121 1 [23lT1) that for one pair of simple turning points the momentum spectrum of 
the pair creation rate (I3T1) in a semiclassical saddlepoint approximation is equal to 


fermion 

/7 k 
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For the bosonic case we can find 


scalar 

"k 



(45) 
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by using the substitution discussed above. We thus find that analogous to the case of two-component 
fields in flat space-time the difference between fermions and bosons is a factor of the form cxp(K[(ij p )). 


4 Pair creation rate of a constant electric field in dS spacetime 


In this section we derive the pair creation rate for a constant electric field, which is described in dS 4 by 



(47) 


This is due to the fact that a co-moving observer, with a four-velocity i/ /J , would measure an electric field 
of 


= m v /> = aE6 z p , 


(48) 


which leads to a constant field strength since £ ; , E 1 ' 
For convenience we introduce 
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where A = jp represent the electric field divided by the Hubble rate as usually in cosmological spacetime 
and y — ^ is the mass term divided by the Hubble rate. These two quantities represent the electric and 
gravitational contribution respectively. We can now check in which limit the semiclassical approximation 
holds. It requires the rate of change of the background to be small in the asymptotic future, that is 


Kwf 

,^ 2 (?])} 



and 
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being small, we see that this is the case when ^ » 1. Comparing to the bosonic case we see that in the 
limit ju » 1 the "-2" in (l23l > is negligible, so that the fermionic and bosonic pulsation are the same. To 
compute the pair creation rate we first have to compute the integrals Ko{r/) and K\ (q) defined in (l30l) and 
(l38l > respectively. The value of ij at the turning point ( 1441 is found to be 


7; -—J-• (51) 

The imaginary parts of Ko(rj p ) and K \ ( ij p ) are the only ones contributing to (l45l >. One can show that 


Im[/f 0 (? 7 p )] = -n[u- jA^8(k z A), 

(52) 

Im[/fi(77')] = 0, 

(53) 


lr The detailed intermediate steps can be found in Eqs. (32)-(38) of (2D. Observe that in the current paper the integration contour 
is closed in the lower imaginary half plane because of opposite convention for the phases in 63- 












Where 0(x) is the Heaviside step function. It was introduced since the real part of the turning point (BTT) has 
to be negative for the turning point to be inside of the closed contour which is needed in the approximation 
of (l43l) . We thus find that only pairs with a momentum in '-direction which has the same sign as A are 
produced in the semiclassical limit. This is what was called pair production in “screening” direction or 
“downward” tunneling in (9). The case of opposite sign, called “upward” tunneling is suppressed in the 
semiclassical limit. 

K\ (r]p) was the only difference between bosons and fermions and is not contributing to the number of pairs 
created. Thus we find that the number of pairs in the semiclassical limit for both bosons and fermions is 
given by 


«k = exp 




9(k z A). 


(54) 


We will close this section by performing the flat spacetime limit and making the relation with the bosonic 
case of S3 explicit. The definition of the pair production rate is 


T = 


(2tt) 3 V 


f 


dk/ik. 


(55) 


where V = a(r/) 4 dri is the unit four volume of the spacetime. As in ED, an estimate for the moment when 
most of the particles are created can be found by analyzing when the adiabaticity is violated. This gives 

0 -T- (56) 
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Using this the /^-integral can be changed into a time integral. Going to spherical coordinates the k- integral 
can be performed. Putting everything together, one finds 
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We can compare this with Eq. (2.33) of El, that we reproduce here 
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We see that in the semiclassical limit /u » 1 and A » 1 the expressions are equal. As in the bosonic case 
the physical number density n of produced pairs at the time rj is given by 


" = ^Xl‘'™ (T)4r = 57?- (59) 

The fact that it is constant shows that the dilution from the expansion of the universe is exactly compen¬ 
sated by the particles created from Schwinger and gravitational particle creation. Hence in the semiclas¬ 
sical limit, the population of fermions is always dominated by the particles created within a Hubble time. 
The vacuum decay rate is defined for fermions as 

Y = log(l - |/? k | 2 ). (60) 


The limit in which the Hubble parameter is negligible compared to the gravitational and electrical strength 
corresponds to the limit to flat spacetime. After some calculations which are analogous to the ones in [f8l 
one can find the Minkowski limit by taking H —> 0, which gives 


(eE) 2 

lim T = —— 7 exp 
h-> o (2 tt) 3 f 


lim Y = V - 

h-» o 4-i l 2 


1 (eE) 2 


(2tt) 3 


nm 

'W\ 

exp 


\eE\ 


( 61 ) 


These are the familiar results for Schwinger pair production in Minkowski spacetime 0. 












5 Conclusions 


In this work, we investigated the fermionic pair creation rate by the combination of an electric and a 
gravitational field in dS 4 . We first presented the basic equations of our setup and derived the corresponding 
Dirac equation in section [2] In section Q] we proposed a semiclassical approximation to compute the 
number of pairs produced. This approximation is based on a saddle point approximation of the integrals in 
©. To make the comparison to the analog bosonic case easier, we first reviewed its computation within 
our formalism. Then we presented the computation of the number of fermions created for a constant 
electric field. It is the main result of this paper shown in (l54l) . The pair creation rate, in this limit, is the 
same in the fermionic case as in the bosonic case. This equivalence between fermions and bosons in the 
semiclassical limit occurs also for one-component fields in flat spacetime when there is only one pair of 
turning points (see e.g. 1201 ). The limit to flat Minkowski spacetime is presented in ( 16 1 1 ) and agrees with 
the usual expression for the Schwinger effect. 

Going beyond the semiclassical limit is possible but is out of the scope of this paper. To do this one 
can compute a more general quantity: the fermionic induced current. The current is proportional to the 
number of pairs created in the semiclassical limit but allows to explore a regime where the notion of 
adiabatic vacuum and of particles does not necessarily exist. It is a more precise way of describing the 
Schwinger effect because its definition does not depend on the time of creation of the pairs and the rough 
estimate (l56l > can be avoided. A calculation of the current induced by an electric field in dS 2 is performed 
in l24l. 

Throughout this paper the electric and gravitational fields have been assumed to be external sources. One 
of the next steps could be to consider backreaction effects of the newly created particles to the external 
electric and gravitational fields. To do so, for the electromagnetic side, the electric current needs to be 
computed and be plugged into the equivalent of the Maxwell equation in curved spacetime. For the 
gravitational side, the density of particle needs to be plugged in the Friedman equation. Another direction 
could be to look for cosmological application e.g. magnetogenesis, baryogenesis and the relation between 
dark energy and dark matter in the context of neutrino physics. 
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